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ABSTRACT 

We used long duration, high quality, unresolved (Sun-as-a star) observations collected by the ground based network BiSON and 
by the instruments GOLF and VIRGO on board the ESA/NASA SOHO satellite to search for solar-cycle-related changes in mode 
characteristics in velocity and continuum intensity for the frequency range between 2.5 mHz < v < 6.8 mHz. Over the ascending 
phase of solar cycle 23 we found a suppression in the p-mode amplitudes both in the velocity and intensity data between 2.5 mHz 
< V < 4.5 mHz with a maximum suppression for frequencies in the range between 2.5 mHz < v < 3.5 mHz. The size of the amplitude 
suppression is 13±2 per cent for the velocity and 9+2 per cent for the intensity observations. Over the range 4.5 mHz < y < 5.5 mHz 
the findings hint within the errors to a null change both in the velocity and intensity amplitudes. At still higher frequencies, in the 
so called High-frequency Interference Peaks (HIPs) between 5.8 mHz < v < 6.8 mHz, we found an enhancement in the velocity 
amplitudes with the maximum 36±7 per cent occurring for 6.3 mHz < v < 6.8 mHz. However, in intensity observations we found a 
rather smaller enhancement of about 5±2 per cent in the same interval. There is evidence that the frequency dependence of solar-cycle 
velocity amplitude changes is consistent with the theory behind the mode conversion of acoustic waves in a non-vertical magnetic 
field, but there are some problems with the intensity data, which may be due to the height in the solar atmosphere at which the VIRGO 
data are taken. 

Key words. Sun: oscillations - Sun: surface magnetism - methods: observational - methods: data analysis 



■ 1 . Introduction 119481: lUlmschneider & NarainI Il990h . Simulations show that 

photospheric oscillations can produce chromospheri c shocks 

>^ The Sun shows a rich spectru m of solar oscillation modes usu- through upward propagation (ICarlsson & SteinI [l997h . but the 

H ; ally referred to as p-mode s Ceibacher & SteinI |1971|; |Ulrich| role this plays in chromospheric heating is still under debate 

P. , 1 19701; iLeighton etaL||1962|; [Deubner||1975|). They are stochas- dFossum & Carlssonll2005l;lKalkofenll2006l) . 
tically excited and intrinsic ally damped by turbulent convection 

dGoldreich & Keeley||1977|) resulting in a significant variation in Through observations and with the help of analyses a va- 

their power with time. In addition, however, there is some evi- riety of types of oscillation modes with frequencies v < v„c 

dence t hat highly energetic eveiits caii excite both high-degree were recendy detected. It was claimed that these are slow- 



modes jKosovichev & Zharkovalll998l) and even global modes magnetoacoustic waves generated by mode transmission dCallvl 
dFoghzzo et alJ|1998|;|Karoff 2003). 1995; Cally & Goosens 2008i) propagating in the atmosphere, 



It was predicted that only those waves with frequencies because the inclined ma gnetic field supposedly r educed the ef- 
below the acoustic cut- off frequency (measured in the Sun fectiv e acoustic cut-off dPe Pontieu et alJ 120041 : iJefferies et al.l 
as Vac ~ 5.3 mHz, e.g. |jimenezll2006|) would be trap ped in- l2006h . The transfer of energy between internal acoustic and ex- 
side the acoustic cavity of the Sun ( Balmfor53ll990l) : while ternal slow magnetoacoustic waves is a function of frequency. 



higher frequency waves, above Vac, could travel freely in thickness of the interaction region (where the Alfven speed is 
the solar atmosphere. Hence th e latter became a possible comparable to the sound speed) and angle between the wave 
source for chromospheric heating dAlfvenll 1 9471 : ISchwarzschildl vector and the magnetic field vector (attack angle). A narrow 



2 



Simoniello et al.: Velocity and intensity amplitude changes 



attack angle, lower frequency and thin interaction region yields 
enhanced transmission from internal acoustic waves to slow ex- 
ternal magnetoacoustic waves (Schunker & Cally 2006). 

This mechanism is the preferred one to explain the lack 
of acoustic power observed in sunspots compared to the 
quiet Sun between 2 mHz < v < 5 mHz both in inten- 
sity and velocity power maps (IWoods & Craml[T98 it iLites et al.l 
[1982: Brown et al. 19921) . Additionally, an acoustic enhance- 
ment above 5 mHz in the velocity amplitudes has been found 
in and/or around sunspots and active regions, mainly in pixels 
with intermediate and weak magnetic field strengths. Nothing 
similar w as found in the continuum intensity (Braun et al. 1992; 
Hindman et al..l998^,Thomas & Stanchfield 2000; Jain & Haber, 



2002ir 



Several questions are still open: Are there any localized 
sources that enhance the acoustic emissivity? Is the same mech- 
anism behind the acoustic absorption and enhancement? Why is 
there a lack of enhancement in continuum intensity data? Are the 
waves incompressible? What, if any, is the influence on the data 
of the height in the solar atmosphere at which the observations 
are localized? In this paper we attempt to address these issues by 
investigating the induced solar cycle changes in velocity and in- 
tensity mode amplitudes for frequencies between 2.5 mHz < v < 
6.8 mHz. We do not limit ourselves to a study of the phenomena 
in sunspots. Instead we decided to use integrated sunlight obser- 
vations, because, although the integr ation over the solar surface 
wiU be dominated by th e quiet Sun, jUnnoll 1 9591: IStenflol[T98l 
ISanchez-Almeidall2004l) mode conversion could take place be- 
cause of the network a nd strong localized mag netic fields with 
strengths of up to 2 kG (^S a nchez- Almeidal2004l) . The simultane- 
ous observation of both absorption and enhancement could be a 
hint that the same physical mechanism is behind them. Equally, 
if we get evidence that at any one time only one of the two phe- 
nomena is observed, this would support the theory of differing 
mechanisms for two processes. 

Furthermore integrated sunlight observations are dominated 
by low-^ modes. The long duration, high quality data available in 
integrated sunlight measurements have already demonstrated the 
interactio n between acoustic waves and the solar magnetic activ- 
ity cy cle telswort h et al.lll993HChaplin et al."2000: Komm et al.l 
200C; Ilimenez-Reves et al.l l2003l ISalabert & Jimenez-RevesI 
20061) . These modes are essentially vertical at the surface and 



so the inclination of the field at the surface is equal to the attack 
angle. 

An additional real novelty of this investigation is the char- 
acterization of solar cycle changes in the velocity and intensity 
mode amplitudes above the acoustic cut-off, and for the first time 
the results from integrated sunlight data are compared with the 
theoretical predictions of mode conversion. Our findings point 
to acoustic suppression in the p-band (2.5 mHz < v < 4.5 mHz), 
and acoustic enhancement in the high frequency band above 
5.8 mHz. 



2. Integrated sunlight observations 

2. 1 . Ground-based network 

The Birmingham Solar Oscillation Network (BiSON) has been 
operating since 1976 and makes unresolved solar disk observa- 
tions. It consists of six resonant scattering spectrometers, which 
perform Doppler velocity measurements in integrated sunlight 
on the K Fraunhofer line at 7699A (Chaplin et al. 1995). The 
data are dominated by the Doppler variations from the low- 
degree i modes. The nominal height in the solar atmosphere at 



which the measurement is made is ^^2 60 km above the photo- 
sphere (e.g. Ijimenez-Reves et aDl2007h as defined by r = 1 at 
500nm . The data provided by the BiSON network for this in- 
vestigation cover the period from April 1996 until July 2006. 

2.2. Satellite instruments 

The Global Oscillation at Low Frequency (GOLF) instrument 
on board the SOlar and Heliospheric Observatory (SOHO) satel- 
lite is devoted to the search of low-degree modes. It works by 
measuring the Doppler shifts of the Na line at 5889A(D1) and 
5896A(D2) jGabrielet al.lll995h . Due to the malfunctioning of 
the polarization system, GOLF is working in single-wing con- 
figuration, which has changed during the 10 years of observa- 
tions as follows: 1) 1996-1998 in the blue-wing configuration; 
2) 1998-2002 in the red-wing c onfiguration; 3) 2002 until today 
in the blue-wing configuration jUlrich et alj|2000l: iGarcfa et al.l 
l2005h . Due to the formation heights of the spectral lines GOLF 
observes at x330 km in the blue-wing configuration, but at 
^ 480 km in the red-wing configuration ( iJimenez-Reves et"al1 
l2007h . For this survey we will use only the first blue-wing period 
together with the red-wing. The data provided by the GOLF in- 
strument for this investigation cover the period from April 1996 
until July 2002. The optimum calibration of the second blue- 
wing period is still under investigation. 

The Variability of Solar IRradiance and Gravity Oscillation 
(VIRGO) instrument is made of three Sun PhotoMeters (SPM) 
at 402nm (blue), 500nm (green) and 862nm (red) that look at 
the Sun as a star (Frohlich et al. 1995)- In this paper we refer 
to the three VIRGO channels as VIRGO(red), VIRGO(green) 
and VIRGO(blue). Each one of the three different signals ob- 
serves the Sun-a s-a-star oscillations a t different heights in the so- 
lar photosphere (iJimenez et al.ll20()5l) . Response functions (RFs) 
are a powerful tool for the analysis of the information content 
and diagnostic potential of spectral lines. These functions mea- 
sure the reaction of the line profile when the atmosphere is lo - 
cally perturbed at a given height (ISocas-Navarro et al ] ll998l) . 
Figure 1 shows the RFs for the three VIRGO channels in their 
respective colors. Although the RFs are broad, the observational 
height can be approximated to be at the maximum value of the 
RFs for the three channels. Therefore we can state that the blue 
channel observes at x!-26km, the green channel at w-lOkm and 
the red channel at x!-i-llkm. Furthermore the blue and green 
channels are more sensitive to the lower part of the solar pho- 
tosphere, whereas the red channel is more sensitive to slightly 
higher photospheric layers. The data provided by the VIRGO 
observations cover the full solar cycle 23. In this survey we use 
the data provided by all three VIRGO channels. The data pro- 
vided by VIRGO for this investigation cover the period starting 
from April 1996 until July 2006. 

3. Data Analysis 

3.1. Mean Fourier transform at low frequency 

We are interested in investigating the frequency dependence of 
the velocity and intensity amplitude with the solar activity cycle. 
Therefore we split the entire frequency range into five intervals 
shown in Table 1 . When one studies the solar oscillations spec- 
trum with helioseismic instruments the convective background 
(granulation, mesogranulation, supergranulation) and the effects 
of active regions crossing through the visible solar disk deter- 
mine the noise back ground and hence Um it the sensitivity of our 
measurements (e.g. Lefebvre et al.ll2068h . There is also the im- 
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Fig. 1. The three response functions for the three VIRGO chan- 
nels (blue, green, red) indicated by their respective colors. 

Table 1. The frequency bands used for the data analysis. 



Interval Name 


Frequency range 
(mHz) 


1st p-band 


2.5 < y < 3.5 


2nd p-band 


3.5 < y <4.5 


3rd p-band 


4.5 < y < 5.5 


1st HIPs band 


5.8 < y < 6.3 


2nd HIPs band 


6.3 < y < 6.8 



pact of time gaps in the dataset. We next discuss our strategy 
for limiting their effect. Our main concern was to accurately de- 
termine the yearly averaged velocity amplitude in each of the 
five frequency bands. Instead of taking the data as a single, year- 
long dataset with gaps we chose to work with a large number of 
subsets, each of which has very high data fill of over 70 per cent. 
This reduces artefacts and improves the accuracy of the mean ve- 
locity amplitude. The choice of the subset length was very care- 
fully considered. At low frequency we found that a subset of the 
length of a half a day was the best compromise between improv- 
ing accuracy and maintaining sufficient r esolution in the acoustic 
power spectrum jSimoniello et al.ll2004h . Hence, to track the so- 
lar cycle changes in the velocity or intensity amplitudes within 
the /9-band, we split the timeseries for each year into subsets with 
a length of half a day. We performed a Fast Fourier Transform 
(FFT) for every half day subset and then we integrated the ve- 
locity amplitudes in each of the five frequency intervals. Finally 
we averaged the 730 subsets to get the yearly mean velocity am- 
plitude ([simoniello et al. 2009a). The errors associated to each 
yearly mean velocity amplitudes were determined by taking the 
yearly standard deviation. 

3.2. Mean Fourier transform at high frequency 

The strongest signals detected in integrated sunlight observa- 
tions are from standing waves whose amplitude results from 
the coherent addition of waves in the solar cavity. Traveling 
waves, however, seldom produce a coherent signal in time, 
and therefore we would not expect to find any such signal 
in the power spectrum of integrated sunlight. Despite this, a 
pattern of equally spaced peaks was found in GOLF, BiSON 
and VIRGO above the acoustic cut-off tGarcfa et al.l Il998bt 
IChaplin et al.l 120031 : iJimenez et al.l l2005h . They were named 
High-frequency Interference Peaks (HIPs), and it is suggested 



that this high-frequency spectrum is the result of geometric in- 
terference between traveling waves (lKumadll99 Ot iKumar & Lul 
[1991; Garcia etal. 1998b). These waves have a much lower q- 
factor than is observed for the peaks around 3 mHz. It has been 
shown that the best subset length to observe the HIPs is four 
days. Hence we analyzed the data by splitting a yearly time se- 
ries into subsets of four days' length. We performed the FFT for 
each subset and then averaged them. Figure |2] shows the high 
frequency part of the power spectrum from the GOLF blue and 
red-wing configuration and BiSON, where the color with which 
the GOLF data are shown is indicative of the observing wing. 
The three different signals have diff erent amplitudes du e to the 
different observational heights dSimoniello et al.lbOOSh . In this 
spectrum we can distinguish three different regions: 

- between 5.0 mHz < v < 5.5 mHz unresolved high-« 
p-mode envelopes (odd and even degree) separated by 60 to 
80 juHz depending on the frequency interval; 

- between 5.3 mHz < v < 5.8 mHz mainly no signal for 
GOLF blue wing configuration, few peaks for BiSON and GOLF 
red-wing configuration; 

- above 5.8 mHz in the GOLF blue-wing configuration up to 
ft! 7 mHz a HIPs pattern that consists of equally spaced peaks of 
almost 70 fiHz. In the GOLF red-wing configuration and BiSON 
we can spot few peaks. 

It has been shown that at higher frequency the noise in 
GOLF observations has an in strumental origin ( iGarcfa et al.l 
l2005l;lTurck -Chieze et al I200I and is dominated by photon shot 
noise. The red-wing configuration has a low photon count com- 
pared with that in the blue wing and we expect the noise levels to 
be commensurately higher This may be an explanations for the 
different behavior between the blue and red-wing configurations. 

3.3. Averaged cross-spectrum at high frequency 

We analyzed the GOLF data with the cross-spectrum technique 
to extract the common signal with a better signal-to-noise ratio. 
The averaged cross spectrum (CS) is defined as the product of 
the Fourier transform of one channel multiplied by the complex 
conjugate of the Fourier transform of the second one and then 
averaged over all the sub-series . It enhances any coheren t sig- 
nals between the two data sets (iGarcfa et al ] |l998bl[T999l) . We 
can apply the cross spectrum technique to the GOLF data, be- 
cause the instrument has two photo-multipliers for each channel 
in the red and blue-wing configuration. Therefore we can de- 
termine the cross-spectrum by using two time series from the 
two photo-multipliers in the blue and red-wing configuration. 
Figure |2] shows the cross-spectrum for the blue and red-wing 
configuration: the noise levels have been reduced by a factor of 
better than 10 and the signal-to-noise ratio has improved by a 
factor of 2. The CS will improve the determination of the HIPs' 
amplitude variation with solar cycle. The lack of two simultane- 
ous time-series prevents us from using BiSON data in the HIPs 
band. Therefore we will apply the cross-spect ral analysis only 
to the GOLF observations in the HIPs band dSimoniello et al.l 
I2009bh . 

4. Frequency dependence in velocity and intensity 
amplitudes with solar cycle 

4.1. Soiar cycie changes in p-mode velocity and intensity 
amplitudes 

We analyzed BiSON, VIRGO and GOLF data starting from 
1996 June 25 until 2006 June 25 for BiSON and VIRGO, while 
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Fig. 2. Left-hand panel: BiSON and GOLF power spectrum in the HIPs band. Right-hand panel: GOLF cross-spectrum for the blue 
and red-wing observations in the HIPs band. 



for GOLF we stopped in 2002. As indicated earlier, we are still 
working on the best calibration of the third GOLF velocity seg- 
ment (blue-wing configuration) and do not use the data because 
the absolute value of the power could be changed by the calibra- 
tion. For BiSON and VIRGO we took the year 1996 as reference 
level of activity and will show the results over the full solar cy- 
cle for the three different frequency bands listed in Table 1 . For 
GOLF we made diff'erent choices. Because it is observing at two 
different heights in the atmosphere, for the blue period we de- 
termined the amplitude variation with respect to 1996 and for 
the red-wing configuration with respect to 1998. The sensibil- 
ity to the visible solar disk is also different between both wings 
as it w as theoretically studied ( iGarcia et al.lI998al : lHennev et al.l 
[T999h . 

Figure |3] shows velocity amplitude changes with solar cycle 
from the BiSON (left panel) and GOLF (right panel) observa- 
tions. It shows several interesting features, and Table 2 quantifies 
the maximum value of the solar cycle changes for the different 
frequency bands: 

- BiSON in the first frequency band shows a net decrease of 
about 14+3%. Within the errors, GOLF observations also show 
a suppression of about 13+2% that agrees well with BiSON. 
Furthermore these findings agree well with previous results 
dChaphn et al.1 120001: iKomm et al.l l2000t |Jim6nez-Reves et al.l 
I2OOI 

- in the second frequency band the p-mode amplitude sup- 
pression is sharply reduced compared to the 1st frequency band 
for both observational programs; 

- in the third frequency band we find almost null change 
in the sign of the variation from the GOLF observations. For 
BiSON we cannot provide any estimate of the solar cycle 
changes in this band due to an instrumental artefact that dom- 
inates any solar cycle changes induced in /9-mode velocity am- 
phtudes. 

We also used VIRGO observations to investigate /9-mode 
amplitude variation in intensity observations. Figure |4] shows p- 
mode amplitude variation from intensity observations: 

- in the first frequency band we found a /j-amplitude suppres- 
sion of 9±2% from all the three channels. Their values indicates 
a slight weaker suppression in this frequency band compared to 
the findings in the velocity observations; 

- in the second frequency band the variation is slightly nega- 
tive; 



- in the third frequency band we find a null change. 

The findings obtained from the velocity observations high- 
light a frequency dependence of the size of the suppression over 
the ascending phase of solar cycle 23. Indeed the findings shows 
a maximal suppression in the first p-band, which becomes less 
pronounced in the second and vanishes in the third p-band. In 
intensity observations, we found the maximal suppression in the 
first p-band, but in the second and third the size of the suppres- 
sion is almost null in both p-bands (Table 2). Nevertheless the 
common features between velocity and intensity observations is 
that the maximal suppression is observed in the first p-band and 
an almost null suppression in the third p-band. 

4.2. p-mode correlation with actvity proxy 

Velocity observations seem to suggest a solar cycle dependence 
of the size of the suppression for two out of three with the pro- 
gression through solar cycle. Therefore we decided to analyze 
the yearly amplitude changes as a function of two well-known 
proxies of global surface activity: the International Sunspot 
Number (ISN) and Mgll H and K core-to-wing ratio. We used 
these two activity proxies, because they are sensitive to differ- 
ent aspects of the magnetic flux. The ISN is sensitive only to the 
strong polarized magnetic flux component, while the Mgll core- 
to-wing-ratio is sensitive to the entire solar disk and therefore 
also includes the co ntribution from bright faculae and network 
dViereck et al.ll200n) . Furthermore was shown that the Mgll ac- 
tivity proxy provides the most consistent description of ;>mod e 
frequency shifts over three solar cycles (IChaplin et al ] l2007h . 
Averages of each proxy were computed over the same periods 
over which the data were analyzed to derive each of the /9-mode 
yearly mean velocity amplitudes. Due to the agreement between 
BiSON and GOLF in the first two p-bands and because ppres- 
sion in the VIRGO observations is almost null already in the sec- 
ond /7-band (we will explain the reason in sec. 5.1), we decided 
to calculate the correlation between the GOLF observations and 
the two activity proxies. Figure|5]shows p-mode velocity ampli- 
tudes variation as seen by GOLF red-wing configuration over the 
period 1996-2001 compared to the changes observed in ISN and 
Mgll. It clearly shows that the size of the p-mode velocity am- 
plitude suppression increases with increasing level of magnetic 
activity for both activity proxies in two out of three p-bands,. 
We then calculated the correlation coefficient with the Spearman 
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Table 2. Maximum amplitude of the p-mode power variation as seen from different observational programs. Positive values indicate 
enhancement while negative values indicate suppression to the reference level. 



Maximal Power Variation (%) 



Interval 


BiSON 


GOLF 




VIRGO 




Name 




RED 


RED 


GREEN 


BLUE 


1st p-band 


-14± 3 


-13± 2 


-10± 2 


-9± 2 


-8±2 


2nd p-band 


-4±3 


-7+2 


-1±2 


-1±2 


-1+2 


3rd p-band 


II 


2+2 


1±2 


0±2 


3±2 



-20 
1994 



. : 2.5<i/<3.5 mHz 
; 3.5<l'<4.5 mHz 




Blue wing 



-20 L . . I ^ 

1994 1995 



: 2.5<i'<3.5 mHz 
: 3.5<l/<4,5 mHz 
: 4.5<l/<5.5 mHz 



Red wing 




2000 
Yeors 



Fig. 3. Left-hand panel: solar cycle changes in the /9-bands from the BiSON observations. Right-hand panel: solar cycle changes in 
the /?-bands from the GOLF blue and red-wing observations. 




Fig. 4. Solar cycle changes in the p-bands from the VIRGO red (left), green (center) and blue (right) observations. 



Table 3. Spearmann rank coefficient between the yearly ampli- 
tude variation in the three p-bands and the two activity proxies 
ISN and Mgll. 



Interval 


ISN 




Mi 


;n 


Name 


Corr 


Rk 


Corr 


Rk 


1st p-band 


-80 


20 


-80 


20 


2nd p-band 


-80 


20 


-80 


20 


3rd p-band 


40 


60 


40 


60 



rank formula, because it also provides the corresponding proba- 
bility of a null correlation. Table 3 indicates a coiTelation of 80% 
in two out of the three p-bands for ISN and Mgll activity proxies. 
These findings agree with previous results ( iJimenez-Reves et all 
l2QQ4HChar)hn et al.ll2007h . 



4.3. Solar cycle changes in the HIPs' velocity and intensity 
amplitudes 

We used 11 years of the VIRGO data starting from 1996 June 
25 until 2006 June 25, while for the GOLF observations we 
stopped in 2002. We split the whole time series in subset length 
of four days to observe the HIPs. We performed the averaged 
cross-spectrum analysis by taking two time series in the blue 
and red-wing configuration of the GOLF observations. Then we 
performed the CS in each of the individual subsets and averaged 
the CS over one year. We considered the HIPs frequency band 
between 5.8 mHz < v < 6.8 mHz in two parts each of the width 
0.5 mHz and present the yearly averages of the signal in these 
bands in Table 4. 

Figure |6] shows the HIPs velocity amplitude variation over 
solar cycle 23 for the GOLF and VIRGO observations. We also 
determined the amplitude vaiation for the blue period in the same 
way as for the lower frequency band, because GOLF observes at 
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Fig. 5. Solar cycle changes in the p-bands from the GOLF blue and red-wing observations compared to the ISN (left-hand panel) 
and Mgll (right-hand panel). In green the ISN (left panel) and Mgll (right hand panel) activity proxies. 



Table 4. Maximum amplitude of the HIPs' power variation as 
seen from the GOLF and VIRGO observations. Positive values 
indicate enhancement to the reference level. 



Maximal Power Variation (%) 
Interval GOLF VIRGO 

Name RED RED-GREEN 

1st HIPs band 18±3 3±2 

2nd HIPs band 36±7 5±2 



Table 5. Spearman rank coefficient between the yearly ampli- 
tude variation in the two HIPs bands and the two magnetic ac- 
tivity proxies ISN and Mgll. 



Interval 


ISN 




M 


gll 


Name 


Con- 


Rk 


Corr 


Rk 


1st HIPs band 


so 


20 


I 


20 


2nd HIPs band 


80 


20 


1 


20 



two different heights in the atmosphere, that is, we determined 
the blue period with respect to the data from 1996 and those 
of the red-wing configuration with respect to those of 1998. We 
found an increase of acoustic emissivity of about 18+3 per cent 
over the period 1998-2001 for the first HIPs band, while ^ 36+7 
in the second HIPs band. This finding seems to point to an in- 
crease of acoustic power with increasing frequency. Based on 
these results, we decided to look for solar cycle changes in the 
HIPs band by using the continuum intensity data provided by 
the VIRGO Sun photometers. The signal-to-noise ratios are poor 
and hence we decided to apply the CS technique to the various 
different color channels of the VIRGO observations, in order to 
enhance any coherent signal. We improved the signal-to-noise 
ratio by a factor of 2. 

Figure |6] (right panel) shows the HIPs intensity amplitude 
variation over solar cycle 23 between the blue and red chan- 
nels. As we can see we find a slight indication of an increasing 
acoustic emissivity over the solar cycle within the errors, but the 
size of the enhancement is pretty small. It is important to under- 
line that the CS analysis performed between blue and green and 
green and red shows up with a pretty similar behavior as the blue 
and red channel. 

4.4. HIPs correlation with actvity proxy 

The findings from the GOLF red-wing configuration seems to 
point to an increasing acoustic power with increasing magnetic 
activity cycle, and the size of the effect increases with frequency. 
Although, we have only four points, we decided to analyze the 
yearly mean amplitude changes in the two HIPs bands as a func- 
tion of our chosen proxies of global surface activity: the ISN 
and and Mgll H and K core-to- wing ratio. Figure |7] shows the 



amplitude variation in the two HIPs bands compared to the ISN 
variation over solar cycle 23 in the left panel, while the right 
panel shows Mgll. In both cases it is abundantly clear that the 
acoustic power increases with increasing solar magnetic activ- 
ity. We calculated the Spearman rank coefficient as in Table 5. 
The results point to a high correlation with the ISN and a cor- 
relation of 1 with Mgll. However, we should emphasize that we 
used only four points out of 11 of the full solar cycle for the two 
determination. 



5. Discussion 

5.1. Observational evidences of mode conversion 

We investigated solar cycle changes in acoustic power with fre- 
quency and with the progression through the solar cycle. In the 
p-bands we found an acoustic power suppression that strongly 
decreases with increasing frequency. Furthermore the GOLF and 
BiSON observations show a good agreement in the strength of 
the effect within the errors. In the HIPs bands we found instead 
an acoustic power enhancement with the solar cycle. Therefore 
we decided to assemble these findings in a composite plot as 
FiglS] shows. It clearly visualizes the change in the sign of the 
variation with increasing frequency. We built up the same com- 
posite plot by using intensity observations. It shows that for in- 
creasing frequency the acoustic enhancement observed in inten- 
sity data is pretty small. The same frequency dependence of 
/9-mode velocity amplitudes was found in MDI observations, 
where p-mode a mplitudes in magne tic regions and quiet Sun 
were compared (Ijain & Habeijl2002h . We now consider an ex- 
planation for the variations. 

Mode conversion has been proposed as a possible mecha- 
nism for describing the p-mode acoustic energy absorption in 
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Fig. 6. Left-hand panel: solar cycle changes in the HIPs bands from the GOLF blue and red-wing observations. Right-hand panel: 
solar cycle changes in the HIPs bands from the VIRGO observations. 




Fig. 7. Solar cycle changes in the HIPs bands from the GOLF blue and red-wing observations compared to the ISN (left-hand panel) 
and Mgll (right-hand panel). In green the ISN (left panel) and Mgll (right hand panel) activity proxies. 



sunspots and active regions (ICallvll2003t ICrouch et al.ll200l . A 
strong surface magnetic field allows the initial internal acous- 
tic wave energy to be split into fast and slow magnetoacoustic 
branches near the Zeq equipartition layer where a ^ c. This split- 
ting has to satisfy the following conservation wave energy: fast 
energy + slow energy = total energy. If the magnetic field is 
enough inclined to lower the acoustic cut-off frequency in low fi 
plasmas, the transmitted slow magnetoacoustic wave may prop- 
agate into the solar atmosphere (ramp effect). If the acoustic cut- 
off frequency has not been reduced, the slow magnetoacoustic 
wave is instead refracted backwards before re aching the observa- 
tion h eight. The transmission coefficient, T dSchunker & Callvl 
l2006h . is defined as the amount of energy transmitted from inter- 
nal acoustic-to-slow or slow-to-internal acoustic waves during 
mode transmission. It has been shown that for the same mag- 
netic fiel d strength and inclination , T is smaller for higher fre- 
quencies dSchunker & Callvll2006l fig. 15). This implies that (in 
the interior of the Sun) the transmission from internal acoustic 
waves to slow magnetoacoustic becomes less dominant with in- 
creasing frequency. Figure |9] shows the transmission coefficient, 
r, of a 6 mHz internal acoustic wave which transmits into the 
low-j8 region as a slow magneto acoustic wave in the gen eral- 
ized ray theory approximation bv lSchunker & CaTlvl (l2006h . For 
a 2 kG magnetic field strength and 60° inclination the transmis- 



sion coefficient is below k 20%, while for a 1 kG magnetic field 
strength it is below » 40%. Here the ray has a lower turning 
depth of 35 Mm, which corresponds to an to harmonic degree 
( ^ 300 with an angle t o vertical at the surface of less than 10° 
dBloomfield et al.l l2007l fig. 7). This is not exactly comparable 
to the global modes case, where I is very small, but it is the 
current depth limit for the code as it is set up for local helioseis- 
mology. Extrapolation to more vertical waves at the surface (and 
smaller C) brings the dashed and solid lines in Fig. |9] together. 
Our observations found a maximum suppression of power be- 
tween 2.5 mHz < V < 3.5 mHz, possibly caused by slow mag- 
netoacoustic waves traveling upwards. Therefore it may be the 
case according to the mode conversion theory that the field in- 
clination is of the order of 60° to allow the resonant modes to be 
transmitted predominantly as slow-magnetoacoustic waves trav- 
eling outwards from the photosphere. For higher frequencies be- 
tween 3.5 mHz < V < 4.5 mHz and for the same field inclina- 
tion, the transmission coefficient is slightly smaller in agreement 
with the reduced acoustic power suppression found in the sec- 
ond p-band. In the third p-band the transmission coefficient gets 
even smaller which allows the resonant modes to be split almost 
equally between the slow and the fast component. This agrees 
with the almost null suppression found in the third p-band. It 
fails, however, to explain the enhancement of power in the HIPs 
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Fig. 8. Left-hand panel: solar cycle changes in the p and HIPs bands from the BiSON and GOLF red-wing observations. Right-hand 
panel: solar cycle changes in the p and HIPs bands from the VIRGO observations. 
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Fig. 9. The transmission coefficient, T, jSchunker & Callvl2006l) 

against magnetic field inclination, 9, for a 6 mHz wave for a ray 
with a lower turning depth of 35 Mm. The blue curve is the 1 kG 
case and the black curve is the 2 kG case. The dashed line is 
for the inclination away from the wave vector and the solid line 
is for the inclination in the same direction as the wave vector. 
When T is small there is strong conversion of internal acoustic 
waves across the equipartition layer Zeq to fast magnetoacoustic 
waves. 



bands. Above the acoustic cut-off, at the 60° field inclination, 
the internal acoustic waves above 6 mHz retain their fast nature. 
In this case, the two restoring forces acting to generate the fast 
magnetoacoustic mode (gas pressure and magnetic pressure fluc- 
tuations) are in phase and add constructively for the fast mode 
increasing the acoustic power. Therefore, we observe a change in 
the sign of the integrated acoustic power above the acoustic cut- 
off with a maximum occurring between 6.3 mHz< v <6.8 mHz. 

The weak point is the lack of enhancement in intensity ob- 
servations. The possible explanation of the observed difference 
between intensity and velocity observations could be found in 
the position of the equipartition layer {Zeq)- The HIPs are trav- 
eling waves in nature, and therefore in order for mode trans- 
mission/conversion to occur the Zeq layer has to lie above the 
excitation so urces that are located around 160 km below the 
photosphere (iKumar & BasullT999l) . For a 1 kG magnetic field 
strength the position of the Zeq is around 80 km above the pho- 
tosphere. Continuum intensity observations, as we have shown, 
are performed below this point, and therefore mode transmis- 



sion/conversion at the VIRGO observational height has not yet 
occurred. This furthermore explains also why in the second and 
third /5-band of the intensity observations we found an almost 
null size of suppression. 

5.2. Leaking fast and slow MHD waves to magnetic features 

We showed that the strength of the suppression in /9-mode ve- 
locity amplitudes is highest when the solar activity is high- 
est, which also holds for the strength of the acoustic enhance- 
ment in the HIPs bands, which is highest when the solar activ- 
ity is highest. Therefore, both p-mode power suppression and 
HIPs power enhancement are correlated with the level of solar 
activity. We have learnt that mode conversion occurs in mag- 
netic field strengths greater than ^ 1 kG and that at a high 
(> 30°) field inclination the transmission coefficient from inter- 
nal acoustic to slow magnetoacoustic becomes much less signif- 
icant. This is extremely important to keep in mind if we want to 
link slow and fast MHD waves to the different solar magnetized 
areas. For example, magnetic regions have already been found 
to be strong absorbers of acoustic power in the frequency inter- 
val 2 mHz < V < 4 mHz, a nd the suppression increases with 
the field strength dJain & Ha ber 2002). Instead high-frequency 
waves seem to show up with a different behavior according to the 
magnetic field strength: the power is still suppressed in regions 
with strong magnetic fields, while areas surrounding them are 
found to be strong enh a ncers of acous t ic emissivity above 5 mH z 
(iHindman et al.1 II998I: iJain & HabeJ 120021: iNagashimal l2007h . 
Why are regions of strong enhancement found around sunspots? 
Intense magnetic fields {^2 kG) with small field inclinations 
are generally confined to the umbra, and this physical condition 
meets the conditions for a larger transmission coefficient. The 
resonant modes (mainly between 2 mHz < v < 4 mHz) are the 
ones that could undergo slow-mode transmission, because the 
Zeq lies below their upper turning point, while in the same mag- 
netic field strength the Zeq for the traveling waves is below their 
excitation sources. As we go away from the umbra, the magnetic 
field strength decreases and a s a consequence the magne tic field 
inclination becomes steeper dSolanki & Schmid3lI993h . These 
two conditions cause a decrease of the transmission coefficient, 
because by decreasing the magnetic field strength the Zeq moves 
above their upper turning points. The peak transmission coeffi- 
cient was actually found to occur close to a magnetic field in- 
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clination of 30° for all frequencies (ICrouch & Callvll2003h . But 
at the same time, the decreasing magnetic field strength moves 
the Zeq above the excitation sources of traveling waves and the 
increasing magnetic field inclination enhances the conversion 
from acoustic to fast-MHD waves. As a result the fast-MHD 
waves can be localized in the outer region of the sunspots, al- 
though recently observations have shown evidences of oscilla- 
tions of 3 minutes inside sunspots (jCenteno et al. 2009). Slow- 
MHD waves are mainly inside the sunspots. This mechanism can 
of course also explain the distribution of o scillations in and/or 
around active regions dJain &Habedl2002l) . Recently, oscilla- 
tions with frequencies between 2 m Hz < v < 8 mHz were found 
concentrated in bright plage areas (iDe Wiin & Mclntoshll2009t 
IO'Shea&Do"vi3l2009h . 

Now that we have found a possible link between the solar 
magnetic features and slow/fast MHD waves, the results from in- 
tegrated sunlight observations can help us to quantify how much 
of the solar surface is affected by mode conversion/transmission. 
Indeed, the strength of the low-frequency acoustic absorption 
and high-frequency acoustic enhancement is a direct measure 
of the total contribution coming from different magnetized ar- 
eas of the solar surface. It has been shown that the solar sur- 
face coverage ar ea of sunspots and plage s over solar cycle 23 
has been of 1% jS anchez- Almeidal l2004l) . Therefore it will be 
extremely important to look for the missing contributi on in the 
network and the intranetwork (Quiet Sun) areas ( Vecchio et al.l 
I2OO9I) . The understanding of the origin of low-frequency acous- 
tic absorption and high-frequency acoustic enhancement is of 
particular interest nowadays, due to the possibility to use the 
induced solar cycle changes on p-mode parameters as precur- 
sor of solar activity. Acoustic waves seem indeed to undergo 
substantial changes befo re the appearance of macroscopic struc- 
ture at the solar sur face ( IChaphn et al. l l2007tlHowee"nni2009t 
ISalabert et al.l2009h . Therefore to correctly use the p-mode solar 
cycle changes as a precursor of activity it is extremely important 
to understand their origin and link them to the specific magne- 
tized areas. 

6. Conclusions 

Low-frequency acoustic absorption in sunspots is very well 
known in the literature and mode conversion is nowadays be- 
lieved to be the likely mechanism behind it. Localized sources 
of excitation have been proposed as a further mechanism to 
explain the observed high-frequency acoustic enhancement in 
and/or around sunspots. We decided to address this point by us- 
ing 1 1 years of velocity and intensity integrated sunlight obser- 
vations provided by BiSON and VIRGO and 6 years provided by 
the GOLF instrument, in an effort to search for changes to mode 
velocity amplitudes over the ascending phase of solar cycle 23. 
The results are dominated by the magnetism of the quiet Sun and 
only if low-frequency acoustic absorption and high-frequency 
acoustic enhancement occur in a significant area of the solar sur- 
face, we can expect to see evidence for both phenomen. 

In agreement with previous findings, our results confirm 
a decrease of about 13% for resonant modes in the range of 
2.5 mHz < V < 3.5 mHz. With increasing frequency (3.5 mHz 
< V < 6.8 mHz) there is a strong indication of a change of 
sign in the velocity amplitude suppression. We found in veloc- 
ity observations at high-frequency (6.3 mHz < v < 6.8 mHz) an 
acoustic enhancement of about 36+7 per cent, while in intensity 
it was about 5+2. This is the first time that integrated sunlight 
measurements have been characterized by the frequency depen- 
dence of velocity and intensity amplitudes on the solar cycle in 



the HIPs band. Therefore the novelty of this investigation is the 
simultaneous observation of low-frequency acoustic absorption 
and high-frequency enhancement. Because of this finding, we 
can conclude that the same physical mechanism is behind in the 
Quiet Sun. The frequency dependence of the solar cycle changes 
in velocity amplitudes allowed us to compare the observations 
with the theoretical predictions of mode conversion. This com- 
parison seems to point to mode conversion as the likely candi- 
date to explain our findings. This statement could be strength- 
ened by performing intensity observations at different heights in 
the atmosphere, in order to check upon the nature of the waves. 
If we do not observe any high-frequency acoustic enhancement 
with height, this will im ply that they are incompressible waves 
(ICallv & Goosensll2008l) . and we can rule out mode conversion 
as the physical mechanism behind low-frequency acoustic ab- 
sorption and high-frequency acoustic enhancement. The wave 
behavior with height is even more extremely important to in- 
vestigate another interesting issue. While fast magnetoacous- 
tic waves are refracted backwards, slow-magnetoacoustic waves 
can travel in the atmosphere. Thus the latter could play a major 
role in the chromospheric heating (lefferies et al. 2006). 

This opens up a new interesting scenario: these observations 
have shown that the magnetism of the Quiet Sun affects the na- 
ture of the waves with increasing activity, which could even im- 
ply that the chromospheric heating varies with the solar cycle. 
Nowadays, the nature of the Quiet Sun magnetism and its link 
with solar activity as the dependence of chromospheri c heating 
on the state of the s olar activity cycle is hotly debated dCattaned 
1 19991: ISanchez-Alm eida 200^. To work out this it will be ex- 
tremely important to check upon the different contributions com- 
ing from regions which have a strong relatively uniform mag- 
netic field, e.g. sunspots and active regions, with that (if any) 
from regions in which the field direction is very variable, like 
the intranetwork. The results of this further work will help us 
settle the debate about the origin of the quiet Sun magnetism. 
Furthermore the investigation will also be important to check 
the possibility of using the solar cycle changes induced in the 
/7-mode parameters as precursors of solar activity. We hope to 
extend this analysis during the peculiar extended minimum of 
activity cycle 23, in which some frequ ency shifts are not fol- 
lowing the tradition al activity indexes (iBroomhall et al.l [20091: 
ISalabert et aLll2009l) . 

Additionally, the interpretation of the observational findings 
in terms of mode conversion helps us to understand why low- 
frequency acoustic enhancement is essentially located around 
sunspots and active regions. We have shown that the magnetic 
field strength and inclination control the locations where mode 
transmission and/or conversion might occur Transmission from 
acoustic to slow MHD waves between 2 mHz < v < 4 mHz is 
more likely to occur within regions of a strong magnetic field and 
small field inclination. These conditions are mainly satisfied in 
the umbra area. Instead conversion from acoustic to fast MHD 
waves above 5.8 mHz is more likely to occur around sunspots 
and active regions, because it requires a steep field inclination 
and a magnetic field strength of the order of 1 kG. These two 
conditions are satisfied in the sunspot penumbra. The analysis 
carried out in this investigation shows that integrated sunlight 
measurements can investigate the solar-cycle changes induced 
in the p and HIPs band. This has a further implication on astero- 
seismology. The technique developed here can be used to extend 
our understanding of stellar activity (lKaroffll2007l: iKaroff et alJ 
2009i). Recent missions like Kepler aim to track variable stars 
for several years, which will give us the chance to follow their 
magnetic activity cycle by using intensity observations at low 
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and high frequency. The analysis of the /j-mode amplitude be- 
havior with stellar activity cycle will give us a measure of the 
strength of the star activity compared to that of the Sun. 
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